3.214 \(\int \sec ^2(c+b x) \sin (a+b x) \, dx\)

Optimal. Leaf size=34 \[ \frac {\sin (a-c) \tanh ^{-1}(\sin (b x+c))}{b}+\frac {\cos (a-c) \sec (b x+c)}{b} \]

[Out]

cos(a-c)*sec(b*x+c)/b+arctanh(sin(b*x+c))*sin(a-c)/b

________________________________________________________________________________________

Rubi [A]  time = 0.03, antiderivative size = 34, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, integrand size = 15, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.267, Rules used = {4580, 2606, 8, 3770} \[ \frac {\sin (a-c) \tanh ^{-1}(\sin (b x+c))}{b}+\frac {\cos (a-c) \sec (b x+c)}{b} \]

Antiderivative was successfully verified.

[In]

Int[Sec[c + b*x]^2*Sin[a + b*x],x]

[Out]

(Cos[a - c]*Sec[c + b*x])/b + (ArcTanh[Sin[c + b*x]]*Sin[a - c])/b

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 2606

Int[((a_.)*sec[(e_.) + (f_.)*(x_)])^(m_.)*((b_.)*tan[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Dist[a/f, Subst[
Int[(a*x)^(m - 1)*(-1 + x^2)^((n - 1)/2), x], x, Sec[e + f*x]], x] /; FreeQ[{a, e, f, m}, x] && IntegerQ[(n -
1)/2] &&  !(IntegerQ[m/2] && LtQ[0, m, n + 1])

Rule 3770

Int[csc[(c_.) + (d_.)*(x_)], x_Symbol] :> -Simp[ArcTanh[Cos[c + d*x]]/d, x] /; FreeQ[{c, d}, x]

Rule 4580

Int[Sec[w_]^(n_.)*Sin[v_], x_Symbol] :> Dist[Cos[v - w], Int[Tan[w]*Sec[w]^(n - 1), x], x] + Dist[Sin[v - w],
Int[Sec[w]^(n - 1), x], x] /; GtQ[n, 0] && FreeQ[v - w, x] && NeQ[w, v]

Rubi steps

\begin {align*} \int \sec ^2(c+b x) \sin (a+b x) \, dx &=\cos (a-c) \int \sec (c+b x) \tan (c+b x) \, dx+\sin (a-c) \int \sec (c+b x) \, dx\\ &=\frac {\tanh ^{-1}(\sin (c+b x)) \sin (a-c)}{b}+\frac {\cos (a-c) \operatorname {Subst}(\int 1 \, dx,x,\sec (c+b x))}{b}\\ &=\frac {\cos (a-c) \sec (c+b x)}{b}+\frac {\tanh ^{-1}(\sin (c+b x)) \sin (a-c)}{b}\\ \end {align*}

________________________________________________________________________________________

Mathematica [C]  time = 0.10, size = 88, normalized size = 2.59 \[ \frac {\cos (a-c) \sec (b x+c)}{b}-\frac {2 i \sin (a-c) \tan ^{-1}\left (\frac {(\sin (c)+i \cos (c)) \left (\sin (c) \cos \left (\frac {b x}{2}\right )+\cos (c) \sin \left (\frac {b x}{2}\right )\right )}{\cos (c) \cos \left (\frac {b x}{2}\right )-i \sin (c) \cos \left (\frac {b x}{2}\right )}\right )}{b} \]

Antiderivative was successfully verified.

[In]

Integrate[Sec[c + b*x]^2*Sin[a + b*x],x]

[Out]

(Cos[a - c]*Sec[c + b*x])/b - ((2*I)*ArcTan[((I*Cos[c] + Sin[c])*(Cos[(b*x)/2]*Sin[c] + Cos[c]*Sin[(b*x)/2]))/
(Cos[c]*Cos[(b*x)/2] - I*Cos[(b*x)/2]*Sin[c])]*Sin[a - c])/b

________________________________________________________________________________________

fricas [B]  time = 0.48, size = 69, normalized size = 2.03 \[ -\frac {\cos \left (b x + c\right ) \log \left (\sin \left (b x + c\right ) + 1\right ) \sin \left (-a + c\right ) - \cos \left (b x + c\right ) \log \left (-\sin \left (b x + c\right ) + 1\right ) \sin \left (-a + c\right ) - 2 \, \cos \left (-a + c\right )}{2 \, b \cos \left (b x + c\right )} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(b*x+c)^2*sin(b*x+a),x, algorithm="fricas")

[Out]

-1/2*(cos(b*x + c)*log(sin(b*x + c) + 1)*sin(-a + c) - cos(b*x + c)*log(-sin(b*x + c) + 1)*sin(-a + c) - 2*cos
(-a + c))/(b*cos(b*x + c))

________________________________________________________________________________________

giac [F(-2)]  time = 0.00, size = 0, normalized size = 0.00 \[ \text {Exception raised: NotImplementedError} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(b*x+c)^2*sin(b*x+a),x, algorithm="giac")

[Out]

Exception raised: NotImplementedError >> Unable to parse Giac output: Unable to check sign: (4*pi/x/2)>(-4*pi/
x/2)Unable to check sign: (4*pi/x/2)>(-4*pi/x/2)2/b*((tan(a/2)^2*tan(c/2)-tan(a/2)*tan(c/2)^2+tan(a/2)-tan(c/2
))/(-tan(a/2)^2*tan(c/2)^2-tan(a/2)^2-tan(c/2)^2-1)*ln(abs(tan((b*x+c)/2)-1))+(tan(a/2)^2*tan(c/2)-tan(a/2)*ta
n(c/2)^2+tan(a/2)-tan(c/2))/(tan(a/2)^2*tan(c/2)^2+tan(a/2)^2+tan(c/2)^2+1)*ln(abs(tan((b*x+c)/2)+1))+(tan(a/2
)^2*tan(c/2)^2-tan(a/2)^2+4*tan(a/2)*tan(c/2)-tan(c/2)^2+1)/(-tan(a/2)^2*tan(c/2)^2-tan(a/2)^2-tan(c/2)^2-1)/(
tan((b*x+c)/2)^2-1))

________________________________________________________________________________________

maple [B]  time = 1.97, size = 888, normalized size = 26.12 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(sec(b*x+c)^2*sin(b*x+a),x)

[Out]

-8/b/(-4*cos(a)^2*cos(c)^2-4*cos(a)^2*sin(c)^2-4*cos(c)^2*sin(a)^2-4*sin(a)^2*sin(c)^2)/(cos(a)*cos(c)*tan(1/2
*b*x+1/2*a)^2+sin(a)*sin(c)*tan(1/2*b*x+1/2*a)^2+2*tan(1/2*b*x+1/2*a)*cos(a)*sin(c)-2*tan(1/2*b*x+1/2*a)*sin(a
)*cos(c)-cos(a)*cos(c)-sin(a)*sin(c))*tan(1/2*b*x+1/2*a)*cos(a)*sin(c)+8/b/(-4*cos(a)^2*cos(c)^2-4*cos(a)^2*si
n(c)^2-4*cos(c)^2*sin(a)^2-4*sin(a)^2*sin(c)^2)/(cos(a)*cos(c)*tan(1/2*b*x+1/2*a)^2+sin(a)*sin(c)*tan(1/2*b*x+
1/2*a)^2+2*tan(1/2*b*x+1/2*a)*cos(a)*sin(c)-2*tan(1/2*b*x+1/2*a)*sin(a)*cos(c)-cos(a)*cos(c)-sin(a)*sin(c))*ta
n(1/2*b*x+1/2*a)*sin(a)*cos(c)+8/b/(-4*cos(a)^2*cos(c)^2-4*cos(a)^2*sin(c)^2-4*cos(c)^2*sin(a)^2-4*sin(a)^2*si
n(c)^2)/(cos(a)*cos(c)*tan(1/2*b*x+1/2*a)^2+sin(a)*sin(c)*tan(1/2*b*x+1/2*a)^2+2*tan(1/2*b*x+1/2*a)*cos(a)*sin
(c)-2*tan(1/2*b*x+1/2*a)*sin(a)*cos(c)-cos(a)*cos(c)-sin(a)*sin(c))*cos(a)*cos(c)+8/b/(-4*cos(a)^2*cos(c)^2-4*
cos(a)^2*sin(c)^2-4*cos(c)^2*sin(a)^2-4*sin(a)^2*sin(c)^2)/(cos(a)*cos(c)*tan(1/2*b*x+1/2*a)^2+sin(a)*sin(c)*t
an(1/2*b*x+1/2*a)^2+2*tan(1/2*b*x+1/2*a)*cos(a)*sin(c)-2*tan(1/2*b*x+1/2*a)*sin(a)*cos(c)-cos(a)*cos(c)-sin(a)
*sin(c))*sin(a)*sin(c)-8/b/(-4*cos(a)^2*cos(c)^2-4*cos(a)^2*sin(c)^2-4*cos(c)^2*sin(a)^2-4*sin(a)^2*sin(c)^2)/
(-cos(a)^2*cos(c)^2-cos(a)^2*sin(c)^2-cos(c)^2*sin(a)^2-sin(a)^2*sin(c)^2)^(1/2)*arctan(1/2*(2*(cos(a)*cos(c)+
sin(a)*sin(c))*tan(1/2*b*x+1/2*a)+2*cos(a)*sin(c)-2*sin(a)*cos(c))/(-cos(a)^2*cos(c)^2-cos(a)^2*sin(c)^2-cos(c
)^2*sin(a)^2-sin(a)^2*sin(c)^2)^(1/2))*cos(a)*sin(c)+8/b/(-4*cos(a)^2*cos(c)^2-4*cos(a)^2*sin(c)^2-4*cos(c)^2*
sin(a)^2-4*sin(a)^2*sin(c)^2)/(-cos(a)^2*cos(c)^2-cos(a)^2*sin(c)^2-cos(c)^2*sin(a)^2-sin(a)^2*sin(c)^2)^(1/2)
*arctan(1/2*(2*(cos(a)*cos(c)+sin(a)*sin(c))*tan(1/2*b*x+1/2*a)+2*cos(a)*sin(c)-2*sin(a)*cos(c))/(-cos(a)^2*co
s(c)^2-cos(a)^2*sin(c)^2-cos(c)^2*sin(a)^2-sin(a)^2*sin(c)^2)^(1/2))*sin(a)*cos(c)

________________________________________________________________________________________

maxima [B]  time = 0.52, size = 387, normalized size = 11.38 \[ \frac {2 \, {\left (\cos \left (b x + 2 \, a\right ) + \cos \left (b x + 2 \, c\right )\right )} \cos \left (2 \, b x + a + 2 \, c\right ) + 2 \, \cos \left (b x + 2 \, a\right ) \cos \relax (a) + 2 \, \cos \left (b x + 2 \, c\right ) \cos \relax (a) + {\left (\cos \left (2 \, b x + a + 2 \, c\right )^{2} \sin \left (-a + c\right ) + 2 \, \cos \left (2 \, b x + a + 2 \, c\right ) \cos \relax (a) \sin \left (-a + c\right ) + \sin \left (2 \, b x + a + 2 \, c\right )^{2} \sin \left (-a + c\right ) + 2 \, \sin \left (2 \, b x + a + 2 \, c\right ) \sin \relax (a) \sin \left (-a + c\right ) + {\left (\cos \relax (a)^{2} + \sin \relax (a)^{2}\right )} \sin \left (-a + c\right )\right )} \log \left (\frac {\cos \left (b x + 2 \, c\right )^{2} + \cos \relax (c)^{2} - 2 \, \cos \relax (c) \sin \left (b x + 2 \, c\right ) + \sin \left (b x + 2 \, c\right )^{2} + 2 \, \cos \left (b x + 2 \, c\right ) \sin \relax (c) + \sin \relax (c)^{2}}{\cos \left (b x + 2 \, c\right )^{2} + \cos \relax (c)^{2} + 2 \, \cos \relax (c) \sin \left (b x + 2 \, c\right ) + \sin \left (b x + 2 \, c\right )^{2} - 2 \, \cos \left (b x + 2 \, c\right ) \sin \relax (c) + \sin \relax (c)^{2}}\right ) + 2 \, {\left (\sin \left (b x + 2 \, a\right ) + \sin \left (b x + 2 \, c\right )\right )} \sin \left (2 \, b x + a + 2 \, c\right ) + 2 \, \sin \left (b x + 2 \, a\right ) \sin \relax (a) + 2 \, \sin \left (b x + 2 \, c\right ) \sin \relax (a)}{2 \, {\left (b \cos \left (2 \, b x + a + 2 \, c\right )^{2} + 2 \, b \cos \left (2 \, b x + a + 2 \, c\right ) \cos \relax (a) + b \sin \left (2 \, b x + a + 2 \, c\right )^{2} + 2 \, b \sin \left (2 \, b x + a + 2 \, c\right ) \sin \relax (a) + {\left (\cos \relax (a)^{2} + \sin \relax (a)^{2}\right )} b\right )}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(b*x+c)^2*sin(b*x+a),x, algorithm="maxima")

[Out]

1/2*(2*(cos(b*x + 2*a) + cos(b*x + 2*c))*cos(2*b*x + a + 2*c) + 2*cos(b*x + 2*a)*cos(a) + 2*cos(b*x + 2*c)*cos
(a) + (cos(2*b*x + a + 2*c)^2*sin(-a + c) + 2*cos(2*b*x + a + 2*c)*cos(a)*sin(-a + c) + sin(2*b*x + a + 2*c)^2
*sin(-a + c) + 2*sin(2*b*x + a + 2*c)*sin(a)*sin(-a + c) + (cos(a)^2 + sin(a)^2)*sin(-a + c))*log((cos(b*x + 2
*c)^2 + cos(c)^2 - 2*cos(c)*sin(b*x + 2*c) + sin(b*x + 2*c)^2 + 2*cos(b*x + 2*c)*sin(c) + sin(c)^2)/(cos(b*x +
2*c)^2 + cos(c)^2 + 2*cos(c)*sin(b*x + 2*c) + sin(b*x + 2*c)^2 - 2*cos(b*x + 2*c)*sin(c) + sin(c)^2)) + 2*(si
n(b*x + 2*a) + sin(b*x + 2*c))*sin(2*b*x + a + 2*c) + 2*sin(b*x + 2*a)*sin(a) + 2*sin(b*x + 2*c)*sin(a))/(b*co
s(2*b*x + a + 2*c)^2 + 2*b*cos(2*b*x + a + 2*c)*cos(a) + b*sin(2*b*x + a + 2*c)^2 + 2*b*sin(2*b*x + a + 2*c)*s
in(a) + (cos(a)^2 + sin(a)^2)*b)

________________________________________________________________________________________

mupad [B]  time = 5.35, size = 254, normalized size = 7.47 \[ \frac {{\mathrm {e}}^{a\,1{}\mathrm {i}+b\,x\,1{}\mathrm {i}}\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}-c\,2{}\mathrm {i}}+1\right )}{b\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}-c\,2{}\mathrm {i}}+{\mathrm {e}}^{a\,2{}\mathrm {i}+b\,x\,2{}\mathrm {i}}\right )}+\frac {\ln \left ({\mathrm {e}}^{a\,1{}\mathrm {i}}\,{\mathrm {e}}^{b\,x\,1{}\mathrm {i}}\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}\,1{}\mathrm {i}-\mathrm {i}\right )-\frac {{\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}-1\right )\,1{}\mathrm {i}}{\sqrt {-{\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}}}\right )\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}-c\,2{}\mathrm {i}}-1\right )}{2\,b\,\sqrt {-{\mathrm {e}}^{a\,2{}\mathrm {i}-c\,2{}\mathrm {i}}}}-\frac {\ln \left ({\mathrm {e}}^{a\,1{}\mathrm {i}}\,{\mathrm {e}}^{b\,x\,1{}\mathrm {i}}\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}\,1{}\mathrm {i}-\mathrm {i}\right )+\frac {{\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}-1\right )\,1{}\mathrm {i}}{\sqrt {-{\mathrm {e}}^{a\,2{}\mathrm {i}}\,{\mathrm {e}}^{-c\,2{}\mathrm {i}}}}\right )\,\left ({\mathrm {e}}^{a\,2{}\mathrm {i}-c\,2{}\mathrm {i}}-1\right )}{2\,b\,\sqrt {-{\mathrm {e}}^{a\,2{}\mathrm {i}-c\,2{}\mathrm {i}}}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(sin(a + b*x)/cos(c + b*x)^2,x)

[Out]

(exp(a*1i + b*x*1i)*(exp(a*2i - c*2i) + 1))/(b*(exp(a*2i - c*2i) + exp(a*2i + b*x*2i))) + (log(exp(a*1i)*exp(b
*x*1i)*(exp(a*2i)*exp(-c*2i)*1i - 1i) - (exp(a*2i)*exp(-c*2i)*(exp(a*2i)*exp(-c*2i) - 1)*1i)/(-exp(a*2i)*exp(-
c*2i))^(1/2))*(exp(a*2i - c*2i) - 1))/(2*b*(-exp(a*2i - c*2i))^(1/2)) - (log(exp(a*1i)*exp(b*x*1i)*(exp(a*2i)*
exp(-c*2i)*1i - 1i) + (exp(a*2i)*exp(-c*2i)*(exp(a*2i)*exp(-c*2i) - 1)*1i)/(-exp(a*2i)*exp(-c*2i))^(1/2))*(exp
(a*2i - c*2i) - 1))/(2*b*(-exp(a*2i - c*2i))^(1/2))

________________________________________________________________________________________

sympy [B]  time = 165.06, size = 5552, normalized size = 163.29 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(b*x+c)**2*sin(b*x+a),x)

[Out]

Piecewise((x/cos(c)**2, Eq(b, 0)), (-1/(b*sin(b*x)), Eq(c, -pi/2) | Eq(c, pi/2)), (-log(tan(b*x/2) - tan(c/2)/
(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**6*tan(b*x/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*
b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c
/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/
2) - 1))*tan(c/2)**6/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4
*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan
(b*x/2)**2 - b) + 4*log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**5*tan(b*x/2)/(b*tan
(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)
**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + 3*log(tan
(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**4*tan(b*x/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 -
b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(
b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - 3*log(tan(b*x/2) - tan(c/2)/(tan(
c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**4/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b
*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/
2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - 8*log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)
**3*tan(b*x/2)/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b
*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2
)**2 - b) - 3*log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**2*tan(b*x/2)**2/(b*tan(c/
2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4
- b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + 3*log(tan(b*
x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 -
4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan
(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + 4*log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(ta
n(c/2) - 1))*tan(c/2)*tan(b*x/2)/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b
*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*
x/2) + b*tan(b*x/2)**2 - b) + log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(b*x/2)**2/(b*ta
n(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2
)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - log(tan(
b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/
2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 +
4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))
*tan(c/2)**6*tan(b*x/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c
/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) +
b*tan(b*x/2)**2 - b) - log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)**6/(b*tan(c/2)**
6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b
*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - 4*log(tan(b*x/2)
+ tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)**5*tan(b*x/2)/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2
)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2
+ b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - 3*log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1)
- 1/(tan(c/2) + 1))*tan(c/2)**4*tan(b*x/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*t
an(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*ta
n(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + 3*log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(
c/2)**4/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**
2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 -
b) + 8*log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)**3*tan(b*x/2)/(b*tan(c/2)**6*tan(
b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c
/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + 3*log(tan(b*x/2) + tan
(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)**2*tan(b*x/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**
6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b
*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - 3*log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1
/(tan(c/2) + 1))*tan(c/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan
(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2)
+ b*tan(b*x/2)**2 - b) - 4*log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)*tan(b*x/2)/(
b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan
(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - log(
tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(b*x/2)**2/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/
2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2
+ b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) -
1/(tan(c/2) + 1))/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*t
an(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b
*x/2)**2 - b) + 4*tan(c/2)**5/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*ta
n(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2
) + b*tan(b*x/2)**2 - b) + 8*tan(c/2)**4*tan(b*x/2)/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan(c/2)**6 - 4*b*tan(c/2
)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)**2 + b*tan(c/2)**2 +
4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) + 8*tan(c/2)**2*tan(b*x/2)/(b*tan(c/2)**6*tan(b*x/2)**2 - b*tan
(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*tan(b*x/2)
**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b) - 4*tan(c/2)/(b*tan(c/2)**6*tan(b*x/2)**2
- b*tan(c/2)**6 - 4*b*tan(c/2)**5*tan(b*x/2) - b*tan(c/2)**4*tan(b*x/2)**2 + b*tan(c/2)**4 - b*tan(c/2)**2*ta
n(b*x/2)**2 + b*tan(c/2)**2 + 4*b*tan(c/2)*tan(b*x/2) + b*tan(b*x/2)**2 - b), True))*sin(a) + Piecewise((log(t
an(b*x/2))/b, Eq(c, pi/2)), (0, Eq(b, 0)), (log(tan(b*x/2))/b, Eq(c, -pi/2)), (-2*log(tan(b*x/2) - tan(c/2)/(t
an(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**3*tan(b*x/2)**2/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*
tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) + 2*log(tan(b*x/2) - tan(c/2)/(tan(c/2
) - 1) - 1/(tan(c/2) - 1))*tan(c/2)**3/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/
2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) + 8*log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2)
- 1))*tan(c/2)**2*tan(b*x/2)/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*
tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) + 2*log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*ta
n(c/2)*tan(b*x/2)**2/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*
tan(b*x/2) - b*tan(b*x/2)**2 + b) - 2*log(tan(b*x/2) - tan(c/2)/(tan(c/2) - 1) - 1/(tan(c/2) - 1))*tan(c/2)/(b
*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/
2)**2 + b) + 2*log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)**3*tan(b*x/2)**2/(b*tan(c
/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2
+ b) - 2*log(tan(b*x/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)**3/(b*tan(c/2)**4*tan(b*x/2)**2
- b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) - 8*log(tan(b*x
/2) + tan(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)**2*tan(b*x/2)/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(
c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) - 2*log(tan(b*x/2) + tan
(c/2)/(tan(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)*tan(b*x/2)**2/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 -
4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) + 2*log(tan(b*x/2) + tan(c/2)/(ta
n(c/2) + 1) - 1/(tan(c/2) + 1))*tan(c/2)/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*
x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) - 2*tan(c/2)**4/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2
)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) - 4*tan(c/2)**3*tan(b*x/2)/
(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*
x/2)**2 + b) - 4*tan(c/2)*tan(b*x/2)/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c/2)**3*tan(b*x/2)
- 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b) + 2/(b*tan(c/2)**4*tan(b*x/2)**2 - b*tan(c/2)**4 - 4*b*tan(c
/2)**3*tan(b*x/2) - 4*b*tan(c/2)*tan(b*x/2) - b*tan(b*x/2)**2 + b), True))*cos(a)

________________________________________________________________________________________